Spatial coherence of thermal near fields 
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Q^Abstract 

q 

We analyze the spatial coherence of the electromagnetic field emitted by a half-space 
. ^ at temperature T close to the interface. An asymptotic analysis allows to identify 
^^three different contributions to the cross-spectral density tensor in the near-field 

regime. It is shown that the coherence length can be either much larger or much 
, ^^ horter than the wavelength depending on the dominant contribution. 

I PACS numbers: 42.72 (black body radiation); 73.20.M (surface plasmons); 42.25.K 
^ (coherence in wave optics); 07.79.F (scanning near field optical microscopy) 
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. ^ The typical textbook example of an incoher- 
^ent source is a thermal source. From the point 
^ pf view of temporal coherence, its spectrum is 
P^iven by Planck's function and modified by its 
_ ^ emissivity. For usual sources, the emissivity is a 
^ smooth function of frequency. Thus, the spec- 
tral width is usually on the order of the peak 
frequency of Planck's function. From the point 
of view of spatial coherence, a thermal source is 
often assumed to be delta correlated. Yet, an ex- 
act form of the cross-spectral density tensor has 
been derived for a blackbody radiator and it has 



been shown that the spatial coherence length is 
A/2 [1]. These exact results seem to support the 
statement that a thermal source is incoherent. 
Yet , one has to analyze more carefully the prob- 
lem when dealing with a real thermal source. In 
this paper, we consider a source that consists 
of a half space filled with a lossy material at 
temperature T. We are interested in the emit- 
ted field so that we assume that there are no 
other sources. Thus there is no incident radia- 
tion illuminating the sample. Note in particular 
that this is not an equilibrium situation. 



Since we explicitly introduce a model for the 
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source, the emitted field contains evanescent 
waves in the vicinity of the interface. These 
evanescent waves are not taken into account 
when deahng with blackbody radiation. Yet, 
they modify the coherence properties of the 
source in the near field as was shown in [2] . The 
effect is particularly striking if a resonant sur- 
face wave is excited. It has been shown that the 
coherence length can be either much larger than 
the wavelength or much shorter than A/2 close 
to the surface. Temporal coherence is also dra- 
matically modified. For example, the emitted ra- 
diation is almost monochromatic when a surface 
wave is excited [3] . These results were obtained 
using a direct calculation of the field emitted 
by a half-space in the framework of fluctuation 
electrodynamics [4-7]. 

The subject of this paper is to analyze the spa- 
tial coherence of the emitted held by means of 
an asymptotic evaluation of the cross-spectral 
density tensor in the near-field limit (interface- 
detector distance small compared to the wave- 
length). This analysis permits to retrieve the 
properties reported in [2] and yields insight into 
the physical mechanism responsible for these ef- 
fects. We are thus able to identify all the pos- 
sible contributions to the cross-spectral den- 
sity tensor : thermally excited surface plasmons, 
skin-layer currents and small-scale polarization 
fluctuations. We show that to a good approx- 
imation, the sum of these three asymptotic 
contributions coincides with the exact result. 
We obtain different characteristic behaviours 
that vary in accordance with the dominant 
term. Surface waves such as surface plasmon- 
polaritons or surface phonon-polaritons yield 
long-range spatial coherence on a scale of the 
surface wave propagation length which may be 
much larger than the wavelength when aborp- 
tion is small. On the contrary, skin-layer cur- 
rents and small-scale polarization fluctuations 
lead to a much shorter spatial coherence length 
that only depends on the distance to the in- 
terface. A surprising consequence of this prop- 



erty is that the macroscopic theory of radiome- 
try may be extended into the mesoscopic regime 
insofar as emission is concerned. Note however 
that this conclusion is based on the assumption 
of a local medium. The ultimately limiting scale 
is thus given by the electron screening length 
or the electron Fermi wavelength, whatever is 
larger [8,9]. 



1 Overview 

1.1 Radiation emitted by a thermal source 

In this section, we review the source theory ap- 
proach we use for the computation of the ther- 
mal electromagnetic field [4-7]. We focus on the 
radiation in the vacuum close to a source that 
we model as a linear dielectric with dielectric 
function e{r;u!). The frequency dependence will 
not be indicated explicitly in the following since 
we calculate quantities at fixed frequency (or, 
equivalently, at fixed wavelength A = 27rc/a;. 
The source radiates because it contains a fluctu- 
ating polarization fleld P(r). The spectral den- 
sity of this field is characterized by the cross- 
correlation tensor Sp{ri,r2) that, according to 
the fiuctuation-dissipation theorem [4,5,10,11], 
is given by 

5p^(ri,r2) = ^nu,/fc^r_i '^''^(^1-^2) (1) 

The Kronecker 6"^^ and the spatial (5-function in 
this formula follow from the assumption that 
the dielectric function is isotropic and local. We 
have taken the normal-ordered form for the po- 
larization field spectrum since we are ultimately 
interested in the electromagnetic field measured 
by a photodetector (given by normally-ordered 
field operators [1,10]). The electric field E(r) ra- 
diated by the polarization P(r) is now given by 
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the Green function for the source geometry 

E,{r)= /drXG^^(r,r')P,(r') (2) 

V 3 

where V is the volume of the source, %. e., the do- 
main where Im£:(r') is nonzero according to (1). 
All quantities in (2) arc understood as tempo- 
ral Fourier transforms at frequency uj. The co- 
herence function W^^{ri,r2) of the electromag- 
netic field is now obtained as a thermal average 
of (2), using the polarization spectrum (1). One 
obtains [2,3] 

k Y 

X G"'*{ri,r')G^''{r2,r') (3) 



The problem is now to evaluate this expression 
analytically and to obtain an estimate for its de- 
pendence on the separation s = r2 — ri between 
the observation points. 



1 separation S 





distance z 



fluctuating polarisation 



be expected for this source geometry, the elec- 
tric coherence tensor depends on the distances 
Zi,Z2 of the observers and their lateral separa- 
tion S. For simplicity, we put in the following 
zi — Z2 — z. We also normalize the coherence 
tensor W^^ to its value for ri = r2 in the case of 
blackbody radiation 



W, 



bb 



(4) 



where as usual k — cu/c. As outlined in ap- 
pendix A, we thus get the following expression 
for the spatial Fourier transform of the coher- 
ence tensor 



w'\S, z) 



W'^{S,z) 



d^K 
(27r)2 



,iK-S-22;Im7 



W'^{K) (5) 



where K denotes a wave vector parallel to the 
interface and (K, 7) is the vacuum wave vec- 
tor of a plane wave emitted by the source. The 
tensor w^^(K) is given in appendix A, eq.(A.lO). 
The integration over K in (5) also includes wave 
vectors |K| > k, describing evanescent waves 
the source excites in the vicinity of the interface 
(the quantity 7 is then purely imaginary with 
positive imaginary part). 



Fig.l: Model geometry for a planar source. 

To proceed in the calculation, we focus on the 
simplified geometry shown in figure 1: an infi- 
nite half-space with uniform dielectric constant 
e, separated by the plane z — from the empty 
half-space 2; > 0. For this arrangement, the 
Green tensor is explicitly known spatial 
Fourier transform with respect to the lateral 
separation S = {s^, Sy) = R2 — Ri. Details may 
be found in [10,12,13] and in appendix A. As to 



1.2 Near field coherence function 

In this subsection, the typical behaviour of the 
field coherence function is discussed. We identify 
several distance regimes showing a very different 
behaviour of the lateral coherence function. An- 
alytical approximations for the coherence func- 
tion are deferred to the next section. 

In figure 2 is shown the 'energy density' (the 
trace of the coherence tensor at coinciding posi- 
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Fig. 2. Energy density Trw^^{S = 0, z) vs. distance from a metal surface. Dots: numerical evaluation of the 
integral (5), solid lines: sum of the asymptotic approximations discussed in the text. In the far field, the 
numerically computed value is taken. 

(a) : log-log scale for gold at A = 620 nm (e = -9.3 + 1.2 i) and at A = 400 nm (e = -1.1 + 6.5 i). The 

dielectric constants arc extracted from [14]. Dashed line: l/z^ power law dominating the extreme near field; 
dotted line: exponentially decaying contribution of excited surface modes. 

(b) : log- linear scale, showing the exponentially decaying surface plasmon contribution for gold at 
A = 620 nm. For comparison, the case of tungsten at A = 500 nm is shown where no plasmon resonance is 
found {s = 4.4 -M8i). 



tions) above a metal surface in double logarith- 
mic scale. One observes a strong increase with 
respect to the far field energy density when the 
distance z is smaller than the wavelength. For 
moderate distances z < X, the energy density is 
dominated by an exponentially increasing con- 
tribution [cf. fig.2(b)]. This is due to the ex- 
citation of surface plasmon resonances, whose 
contribution is calculated analytically in sub- 
section 2.1. The other curve in fig. 2(b) shows 
the energy density for the case of tungsten with 
Ree > —1 where no surface mode exists and 
no exponential increase is found. For small dis- 
tances z -C A, the energy density follows a 1/^;^ 
power law ('static limit') that is discussed in 
subsection 2.2. The prefactor of this power law 
involves the imaginary part of the electrostatic 
reflection coefficient lm[{e — 1)/ (e + l)]. The sec- 
ond curve in flg.2(a) illustrates the resonantly 
enhanced energy density for a wavelength where 
Re£ ~ —1. The 'static limit' contribution then 
overwhelms that of the plasmon resonance. 

In figure 3, we show the normalized lateral co- 
herence function at chosen distances from the 
interface. In the far field [plot 3(a)], the coher- 



ence length is A/2, and the coherence function 
the same as for the blackbody field {{sin ks)/ks 
behaviour). This is not surprising since at large 
distances z ^ X, only propagating plane waves 
radiated into the vacuum half space contribute 
to the field. 

When the surface plasmon excitation dominates 
the field energy {z < A), the field coherence 
extends over much longer distances [plot 3(b)]. 
This is because of the weak damping of the plas- 
mon modes in this case. We show below (sub- 
section 2.1) that the coherence length is indeed 
given by the plasmon propagation length. The 
figure also shows that the field is strongly po- 
larized perpendicular to the interface, as is the 
surface plasmon mode. 

At distances z <^ X even closer to the source, 
the field coherence length gets shorter again 
[plot 3(c)]. We show below that in this regime, 
the field behaves as if it was quasi-static (sub- 
section 2.2). This leads to a lateral coherence 
length equal to the vertical distance z from the 
interface and hence much shorter than the wave 
length. We thus find the surprising result that 



4 



thermal near fields have no lower limit in their 
coherence length, as long as the dielectric func- 
tion of the source may be taken as local. 

One might finally ask whether the skin depth 
6 [defined in (23)] is relevant for the radiation 
emitted by a metallic source. This question is 
discussed in subsection 2.3 where we show that 
in the regime S <^ z <^ X, a different power law 
(oc governs the energy density (see fig.4(a) 
below). The lateral coherence behaves similar to 
the static regime z S, however, as shown in 
fig. 4(b) below. 

2 Analytical approximations 



2.1 Plasmon contribution 



It is well known that a dielectric-vacuum in- 
terface supports surface plasmon polariton (or 
phonon polariton) modes provided the dielec- 
tric constant satisfies Ree < — 1 [15,16]. These 
surface modes propagate parallel to the surface 
with a wave vector i^Tpi and are exponentially 
localized in the direction perpendicular to the 
interface. In addition, if there are losses in the 
dielectric, the propagation parallel to the inter- 
face is damped which may be described by a 
nonzero imaginary part of i^Tpi. Mathematically, 
we obtain the plasmon dispersion relation by 
identifying the poles of the transmission coeffi- 
cients t^ {n = s, p) as a function of the wave vec- 
tor K. Only the p-polarization (magnetic field 
perpendicular to the plane of incidence) gives a 
pole at the (complex) position 



pi 



1' 



ReKpi, ImKpi > 0(6) 



The plasmon pole shows up as a sharp peak 
when the integrand w'^^(K) in (5) is plotted as 
a function of the lateral wave vector magnitude 
K [cf. eq.(A.lO)]. This suggests that we get a 



good approximation to the plasmon contribu- 
tion by taking slowly varying terms outside the 
integral and evaluating the pole contribution by 
contour integration. For example, the denomi- 
nator of the \tp\^ transmission coefficient may 
be approximated as 



Mel 



1^7 + 72 1 



1 e 



rim 



(7) 



where 7 = ^{K) and 72 are the vertical wave 
vector components above and below the inter- 
face. It is essential for the contour integration 
to work that one expresses the absolute square 
on the left hand side as the imaginary part of 
an analytic function of K (right hand side). 

It is easily checked from (A. 10) that the trace 
of w^^CK) only depends on the magnitude K of 
the lateral wave vector K. The integration over 
the angle between K and S therefore gives 



Ttw'^{S,z) 
^KdK 



Jo{Ks) e-'''^"'^Trw''{K) 



2n 



(8) 



where Jo(") is the ordinary Bessel function and 
s — [S[. The individual tensor components also 
involve Bessel functions J2{Ks), as discussed 
in appendix B.l. The integration over K may 
be done using the identity (B.2) proven in ap- 
pendix B.2.1. The diagonal elements of the co- 
herence tensor finally take the suggestive form 



w"(S,^)ftiCpie-2'^p'"-* 

-Re 
2 



g"[u) 



Ho{u) - H2{u) 



2 



H.{u) + H2{u) + 



4i 
4i 



\e\ReHQ{u) 



37r 


e\ 


2 


2 


£ + l[7/2 



(9) 
(10) 

(11) 
(12) 

(13) 
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skin depth 8= 0.05 X 



0.01 0.02 0.03 0.04 

s/X 

Fig. 3. Normalized lateral coherence functions for three fixed distances z, plotted vs. the lateral separation 
s = |S|. All plots are for a gold surface at A = 620 nm. Dots: numerical evaluation of (5), solid lines: ana- 
lytical approximations discussed in the text. The numerically computed values were used to normalize all 
curves. 

(a) : far field regime z = 10 A. The trace of the coherence tensor is plotted, normalized to its value for s = 0. 
Solid line: free space coherence function sm{ks) / {ks) . 

(b) : plasmon dominated regime z = 0.1 A. The components w^^ and w^^ are plotted, normalized to (the 
numerically computed) w^^{S = 0,2;). 

(c) static regime z = 0.01 A. The components w^^ and w^^ are plotted and normalized as in plot 3(b). The 
solid curve only contains the extreme near-field contribution (17). 

where Kpi — Imji^Kpi) is the perpendicular 

plasmon decay constant and i?o,2(-) = -f^o,^2(') 
are Hankel functions (or Bessel functions of the 
third kind) [17]. The superscripts ||,_L indicate 
the directions parallel and perpendicular to the 
separation vector S in the interface plane. 



Eq. (9) shows that the plasmon resonance gives a 
contribution to the energy density that increases 
exponentially when approaching the source. 
This behaviour is reproduced by the numerical 
evaluation of (A. 10), as shown in fig. 2(b). As a 
function of the lateral distance s, the correlation 
tensor eq.(9) shows damped oscillations whose 
wavelength is fixed by the plasmon wave vector 



Kpi, as shown in fig.3(b). These oscillations can 
be made explicit using the asymptotic form of 
the Hankel function [17] 



|i^pis|»l: 
^"^^P^^^ ^ /^J;'''^''"'""^''""^'^ ^^^^ 

We thus conclude that the propagation distance 
of the plasmon resonance, as contained in the 
imaginary part of K^i, determines the coherence 
length of the field in this regime. For a dielectric 
constant with small imaginary part, the inverse 
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propagation distance is approximately 
{ReKpif Ime 



ImXpi 



P 2(Re£)^ 



< k. (15) 



Thermally excited plasmons thus lead to a spa- 
tially coherent field on a length scale well ex- 
ceeding the vacuum wave length. They also cre- 
ate a net field polarization, as shown by the 
anisotropy of the tensor elements in (10-12) [see 
also fig. 3(b)]. This anisotropy may be under- 
stood from the fact that the coherence between 
points separated by S is created by plasmons 
propagating parallel to this direction, and the 
latter are polarized in the plane spanned by S 
and the normal vector e^- 



2.2 Extreme near field: quasi-static regime 



We now turn to the near field limit z <^ \. 
Inspecting the integrand of (A. 10), one finds 
that in addition to the plasmon resonance, large 
wave vectors K ^ k dominate the integral. This 
is because the exponential cutoff e^"^^^"^"/^^^ ^ 
Q-2zK gg^g effective only ioi K > 1/z ^ k. We 
thus obtain the asymptotic behaviour of the in- 
tegral when we expand the integrand to leading 
order in the limit 1/z > K ^ k. The transmis- 
sion coefficients, e.g., become in this limit 



4| 


72 1 
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\e\ 




X2|£ + l 


2 




(16) 



We perform the integration over K using 
eq.(B.8), as explained in appendix B.2.2 and get 
the following asymptotic form for the diagonal 
elements of the cross-correlation tensor 



w"(S, z) 



a 



xnf 



{kzy 



9\s/z) 



(17) 



g\\u) 
9\u) 



1 - mV2 

;i + kV4)V2 
1 

;i + ^2/4)3/2 

2 - u^/A 

:i +^2/4)5/2 



Cxnf = — Im 



3 Im£ 



32 £ + 1 16|£ + 1|2 



(18) 
(19) 
(20) 
(21) 



The coherence tensor given by (17) shows a 
power law increase 1/z'^ when the interface is 
approached, as plotted in fig. 2(a). It therefore 
takes over compared to the plasmon contribu- 
tion in the 'extreme near field limit' z ^ l//tpi. 
In this regime, the lateral coherence of the field 
is characterized, as shown by the scale func- 
tions g^{s/z) in (18-20), by a lorentzian shape 
whose scale is set by the distance z to the source. 
Hence, the closer one detects the field, the more 
it is spatially incoherent. 

This behaviour may be understood from elec- 
trostatics: in the near field, the electromagnetic 
fields behave as if they were quasi-static because 
they vary on a length scale much smaller than 
the wave length (retardation is negligible). A 
near field detector is thus sensitive to a source 
of the order of nz^, and spatial coherence 
is observed when these areas overlap, hence for 
a separation smaller than the distance z. Simi- 
lar arguments have also been put forward to in- 
terpret subwavelength resolution in optical near 
field microscopy [18,19]. The electrostatic anal- 
ogy may be pushed even further: it is easily 
checked that we get the same result as (17) using 
electrostatic image theory. As a consequence of 
the fiuctuation-dissipation theorem [4,5,10,11], 
we have indeed 



W^^^'(r2,ri) ocImG'^^(r2,ri) 
otIm^Lage(r2;fi,dj), 



(22) 



where G^^ is again the electric Green function. 
The electric field Eimage is created by the image 
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dj of a dipole (polarized along the a;-^-axis) 
at position ri, the image dipole being located 
at the mirror position ri = {xi,yi, —Zi). This 
image dipole field dominates the Green func- 
tion G*^ at sufficiently close distance from the 
source if the electrostatic refiection coefficient 
{e — !)/(£ + 1) has a nonzero imaginary part 
[see eq.(21)]. 

We stress that there is no lower limit to the spa- 
tial coherence length of the near field, provided 
one uses the framework of a local dielectric sus- 
ceptibility. Model calculations for a free elec- 
tron gas confined to a half-space show that this 
framework breaks down at wave vectors K of 
the order of the Fermi wave vector [8,9]. For our 
problem, this corresponds to typical distances of 
the order of 0.1 nm that are difficult to achieve 
for near field probes even in the optical range. 



2.3 Relevance of the skin depth 



It has become clear from the two preceding sub- 
sections that the lateral coherence of near field 
radiation strongly depends on the distance of 
observation to the source. It might have been ex- 
pected that the skin depth shows up in this dis- 
cussion, since it governs the penetration depth 
for electric fields into the metal. We conclude 
our analytical work by identifying the relevance 
of this length scale. 



Recall that the skin depth is given by 

A/27r 



5 = 



Im y/e 



(23) 



For good conductors and low frequencies, the 
dielectric function is dominated by its zero- 
frequency pole 



1(7 



e — 



Equ; 



(24) 



where a is the (possibly frequency dependent) 
conductivity. In particular, one has \e\ ^ 1 
and (5 ^ A in this regime: the source material 
hence approaches a perfect conductor. This im- 
plies that the large K expansion of the trans- 
mission coefficients in (A. 10) has to be recon- 
sidered: while the limit K ^ k may be jus- 
tified, the limit K ^ may be not. We 
find that there exists an intermediate distance 
regime 5 <^ z <^ X (corresponding to wave vec- 
tors k\y/s\ ^ K ^ k) where the coherence 
tensor shows a different behaviour [20]. The ex- 
pansion of the transmission coefficients in this 
regime reads 



4| 


72 


2 


£ 






(25) 



We finally get an isotropic coherence tensor 



^^Xs,.)-^^^/(V.) 



1 



(1 +1^2/4)3/2 



(26) 
(27) 



The skin layer dominated regime is thus char- 
acterized by a l/2;2 power law for the energy 
density. As shown in figure 4(a), the skin depth 
5 separates this regime from the extreme near 
field where a different power law 1/z^ prevails. 
We observe from fig.4(b) and (26) that the lat- 
eral coherence length is equal to the distance 
z from the source, as in the extreme near field 
regime. This is not so surprising since the field 
propagation in the vacuum half space above the 
source is governed by the length scales A and 
z, whatever the smaller, while the skin depth 
is only relevant for the propagation inside the 
source. 
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Fig. 4. Near field coherence in the skin-effect dominated regime -C z <C A, typical for metals in the far 
infrared. Dots: numerical evaluation of the integral (5), lines: analytic asymptotics discussed in the text, 
(a): energy density Trw^^S = 0,z) above a metallic surface with e = —8.26 + ilO^ (gold extrapolated to 
A = 3.7 nm). The solid line is the sum of the asymptotic contributions derived in this paper, (b): normalized 
lateral coherence functions w^^{S, z)/uf^{0, z) (i = x, z) for fixed distance z, plotted vs. the lateral separa- 
tion s = |S|. The X- and z-polarizations differ because the plasmon contribution already comes into play. 
The numerically computed u;^^(0, z) was used to normalize all curves, this is why the analytic correlations 
exceed unity. 



To conclude, we recall that the different contri- 
butions to the correlation tensor originate in dis- 
tinct domains on the 7^-axis in the integral (5). 
The total correlation tensor is therefore given 
by the sum of the surface plasmon, extreme near 
field, and skin-layer contributions. The accuracy 
of this approximation is visible in figs. 2, 4(a). 
Note that in the figures, the numerically com- 
puted far-field energy density has been added to 
get the correct large distance limit. 

3 Conclusion 



In the near field, the spatial coherence of ther- 
mal radiation differs strongly from the black- 
body field. Confined field modes like surface 
plasmon polaritons that propagate along the 
source surface make the field spatially coherent 
over large scales if they dominate the radiation 
density. At close distance (smaller than the skin 
depth), the radiation is dominated by quasi- 
static fields, and the coherence length drops well 
below the wave length, being limited only by the 
(non)locality of the dielectric response of the 
source material. The cross over between these 
regimes is determined by the skin depth and 



the electrostatic reflection coefficient. We con- 
clude that in the near field, macroscopic con- 
cepts like a local emissivity are still meaningful 
at the subwavelength scale, provided coherent 
surface excitations are absent or subdominant. 



The asymptotic forms for the cross spectral 
density tensor obtained in this paper are use- 
ful to characterize thermal noise fields that 
may perturb particles in integrated microtraps 
'mounted' with electromagnetic fields above a 
solid surface [21-28]. The concomitant scatter- 
ing and decoherence of the guided matter waves 
is discussed elsewhere [29,30]. 

The fluctuation electrodynamics used in this pa- 
per enabled us to treat a non-equilibrium situa- 
tion (thermal source in vacuum at T = 0) where 
the fluctuation-dissipation theorem for the elec- 
tric field is not immediately applicable. In par- 
ticular, we neglected the zero-point radiation 
impinging on the interface from the empty half- 
space. The domain of validity of this approxi- 
mation, as well as the calculation of anti-normal 
ordered correlation functions will be the subject 
of future investigations. 
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A Notations for the plane interface 



The Green tensor describing the emission from 
the source z' < into the vacuum half space 
2; > may be written in spatial Fourier space 
as [10,12] 



"'^ ^ / §^2^''''^''~'''^GHK, Z, z') 

(A.l) 



2^072 



M=s,P 



We use bold capitals to denote vectors parallel 
to the interface, e.g., K = [k^, ky, 0). The ver- 
tical components of the wave vectors in vacuum 
and inside the source are, respectively. 



When the Green tensor (A.l) is inserted into 
the integral (3), the spatial integration over R' 
yields a (5-function for the lateral wave vectors. 
The integration over z' is then 



r 1 
J 2Im 



Ini72 



where the convergence is ensured by the posi- 
tive imaginary part of 72. The resulting coher- 
ence tensor is then of the form (5). We use the 
identity 



kHm £ = 2 Im 72 Re 72 



and get after some elementary algebra: 



j, , 37rRe72 
Ak 72 p 



(A.9) 



xEeSe2;|eif¥|t.r (A.IO) 



K2 



lm7 > 



(A.2) 



72 = +Vek'^ - K2 Im72 > (A 3) ^ Components of the coherence tensor 



The polarization vectors for the s- (or TE-) and 
p- (TM-)polarized waves are taken as 



In this appendix, we outline the calculation for 
the components of the coherence tensor. 



_ ^(2) 



K X e^ 



e(2)^ 



Kz--fK 
k 

Kz - 72K 

Vek 



(A.4) 
(A.5) 

(A.6) 



where K is the unit vector parallel to K. Finally, 
with this choice for the polarization vectors, the 
Fresnel transmission coefficients are 



272 
7 + 72' 



272 

£7 + 72 



(A.7) 



B.l Angular integrations 



The only quantities in (A.IO) that depend on 
the angle (/? between the lateral wavevector K 
and the separation S arc the polarization vec- 
tors e^. To simplify the calculation, we choose 
the X-axis parallel to S. We then get the follow- 
ing azimuthal integrals (eq.9.1.18 of [17]) 

J^^^Kscos^ / sm \ ^ ^^^^^^ ^ J2(i^s)(B.l) 
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The integrals with sin cos tp vanish due to par- 
ity. We also note that one also gets nonzero off- 
diagonal elements W^^ , W^^ due to p-polarized 
modes. For simplicity, these are not discussed 
here. 



B.2 Radial integrations 



nentially for |a;| oo there [see eq.(14)]. There- 
fore, closing the integration contour with a half 
circle, the integral is given by the residue at the 
pole X — +q (because Img > 0), and we get 



X dx 
— 



mHi'\qs) (B.5) 



We are left with integrals over the radial wave 
vector K. These are worked out using the defi- 
nitions (A. 6) of the polarization vectors and the 
transmission coefficients (A. 7). 



B.2.1 Plasmon pole 

To find the plasmon contribution, we extract, as 
mentioned in the main text, the pole of the tp 
coefficient and approximate the other factors by 
their values at the pole. The remaining integral 
can be reduced to the following standard form 



xdx 



x^ 



;Jo(xs) 



IB.2) 



for Img > 0, s > 0. To prove this identity, we 
use contour integration. The Bessel function is 
written as ([17], eqs. 9.1.3, 9.1.39) 



Joix) 



1 



ifP(a;)-iJ^^^(e"^^ 



r(l)/ 



x) (B.3) 



This proves (B.2). Taking the imaginary 

part, we find both the trace and the zz- 
component (12) of the correlation tensor (9). 

For the xx- and yy components of the coher- 
ence tensor, we also need the integral (B.2) with 
the Bessel function J2{xs) instead of Jo{xs) [cf. 
eq.(B.l)]. Using the same reasoning as above, 
this integral is transformed into 



1 
2 



X Ax 
^2 _ q2 



Hi'\xs) 



(B.6) 



In addition to the pole at a; = g, we now have 
a contribution from the — 4i/ (,ts)^ singularity of 
the Bessel H2^\xs) function at the origin. This 
singularity lies on the integration path and is 
therefore taken into account by half its (neg- 
ative) residue at x = 0. Combining the latter 
with the residue at a; = g, we get 



/ 



a; da; 



Hi^\xs) ^i7rH^"'{qs) 



x^ 



An 
q2g2 



[B.7) 



where Hq^\x) is the Hankel function. The inte- 
gral may now be written as 

with an integration path running just above the 
negative real axis. The Hankel function is ana- 
lytic in the upper half plane and vanishes expo- 



We may verify the sign of the second term by 
checking that the function (B.7) vanishes in the 
hmit s — > 0, as is the case for the left hand side 
of (B.6). 

B.2. 2 Near field regime 

In the near field regimes K ^ ^lv^| (extreme 
near field) and k\^/e\ ^ K ^ k (skin layer 
dominated regime), the expansions (16, 25) of 
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the transmission coefficients are straightforward 
to obtain. The final integration involves integer 
powers of K times products of Bessel functions 
and exponentials and is performed using the fol- 
lowing identity 

oo 

jiKMKs)^-^-''^^^ (B.8) 

This may be proven starting from the Fourier 
expansion of the Coulomb potential (writing 
k=(K,A;^)) 

and evaluating the integral over the vertical 
wave vector component with contour integra- 
tion (for 2; > 0, a single pole at — i|K| con- 
tributes). The derivatives of (B.8) with respect 
to z and s then provide all necessary integrals. 
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